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Quenching of quantum Hall effect and the role of undoped planes in multilayered 

epitaxial graphene 
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We propose a mechanism for the quenching of the Shubnikov de Haas oscillations and the quantum 
Hall effect observed in epitaxial graphene. Experimental data show that the scattering time of the 
conduction electron is magnetic field dependent and of the order of the cyclotron orbit period, i.e. 
can be much smaller than the zero field scattering time. Our scenario involves the extraordinary 
graphene n = Q Landau level of the uncharged layers that produces a high density of states at 
the Fermi level. We find that the coupling between this n = Q Landau level and the conducting 
states of the doped plane leads to a scattering mechanism having the right magnitude to explain 
the experimental data. 
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Graphene is a two dimensional carbon material which 
takes the form of a planar honeycomb lattice of sp^ 
bonded atoms. Its two dimensional character and the lin- 
ear dispersion relation implies that, close to the charge 
neutrality point, the electrons obey an effective mass- 
less Dirac equation. The properties of electrons in 
graphene, deriving from the Dirac equation, are funda- 
mentally different from those of electrons in standard 
semi-conductors which obey the Schrodinger equation. 
A remarkable example is the quantum Hall effect which 
is quantized with integer plus half values [iL I4] and has 
even been observed at room temperature [3|. Another 
major fact is the large electronic coherence [J, l5|: that is 
observed even at room temperature, and gives the hope 
to produce devices with new properties. 

Graphene can be produced by mechanical exfoliation 
of graphite [l|, I4I which produces either single or few 
layer graphene. Another way of producing graphene is 
by thermal decomposition of the surface of SiC [J, la, la| ■ 
Epitaxial graphene is attractive because it grows on an 
insulating substrate, with a high structural quality and 
very smooth graphene layers [71]. Both methods pro- 
duce graphene samples with spectacular electron coher- 
ence properties. Yet magnetotransport remains puzzling 
in epitaxial graphene. On one hand, at low magnetic 
field, the experimental results indicate electronic mean 
free paths of more than half a micron at 4K [5|, Isj . On the 
other hand Shubnikov de Haas oscillations of the mag- 
netoresistance are weak and the quantum Hall effect is 
not observed. Here we show that the small amplitude of 
the Shubnikov de Haas oscillations is not a signature of 
a low quality graphene. We show, by analysing the ex- 
perimental data, that the results can be explained by a 
scattering time of the conduction electron which is mag- 
netic field dependent, and reduced to the order of the 
cyclotron orbit period. We argue that the conducting 
states in a doped layer can couple to the zeroth Lan- 
dau levels in an undoped layer, which is on top of the 



doped one. As a consequence the conducting electrons of 
the doped layer are subjected to a scattering mechanism 
that increases with magnetic field, because the number of 
states in the zeroth Landau level increases with magnetic 
field. At low magnetic field the scattering time r can be 
long but at stronger field r decreases in such a way that 
LOT ~ 1, where to is the cyclotron frequency. This forbids 
the observation of strong Shubnikov de Haas oscillations. 
We show also that in some limits the magnetoresistance 
increases linearly with the magnetic field. This has been 
observed recently in epitaxial graphene multi-layers and 
the magnitude of the linear magnetoresistance is in quan- 
titative agreement with our model. 

The letter is organized as follows: we present first 
the experimental results, and show that they can be ex- 
plained by a linear dependence of the scattering rate with 
the magnetic field. Based on experimental results we pro- 
pose a simple model of electronic structure and derive the 
equations satisfied by the Green's functions, in the Self- 
Consistent Born Approximation. Finally we analyse the 
electronic structure of this model and show that its mag- 
netotransport properties are in good agreement with the 
experimental observations. 

Summary of previous experimental results 
Electronic properties of epitaxial graphene have been 
analyzed in great detail. They clearly point to the ex- 
istence of a electron doped plane at the SiC/graphene 
interface [9|, |lO| carrying the main part of the current [5|] . 
Since the screening length is of 1-2 interlayer spacing, 
only 2-3 planes are doped, and the other planes are quasi- 
neutral, as shown in infra-red spectroscopy [ill and thus 
comparatively poor conductors. The experimental evi- 
dence that one of the doped planes carries the main part 
of the current suggests that the other doped planes have 
a low conductivity due to the impurities spatial distri- 
bution. The doped and quasi-neutral planes all have the 
characteristics of an isolated graphene plane i. e. massless 
Dirac electrons [5|, l9|, [lO, lUj This was initially surprising 




FIG. 1: Left: Experimental values of resistivity as a func- 
tion of B for a 100/im x lOOO/xm (dashed line), a l/im x 5fj,m 
(dashed dotted line), and a 0.27/im x 6jj,m (fuU line) samples. 
Main panel: longitudinal magnetoresistance pxx', enclosed 
panel: transverse magnetoresistance pxy Right: pxy/pxx 
(which is equal to lot in the one band model) (lower panel) 
and deduced scattering rate (upper panel). Note that the lon- 
gitudinal magnetoresistances pxx have been already published 
in the case of the 0.27pm x Qpm p] and the lOO^jn x 1000pm 
Ql samples. 



since the epitaxial graphene samples studied in [5| and 
ll| consists of several stacked graphene planes. Yet, as 
proven recently, these planes are rotationally stacked [l^ 
and the relative rotation strongly diminishes the effective 
electronic coupling between the plan es as c omp ared to 
the case of AB (Bernal) stacking [ll, [H, [ll El [11]. 
Analysis and interpretation of experimental results 



The magnetotransport coefficients pxx and pxy for 
three samples of various widths (0.27/im, 5/im and 
lOOOfim) are plotted in figure ((T]). Because the quan- 
tum corrections to the conductivity are small [8| a semi- 
classical one band model should be a fair approxima- 
tion. The transverse magnetoresistance pxy increases es- 
sentially linearly with the magnetic field B and is indeed 
consistent with a semi-classical model. Yet, the longi- 
tudinal magnetoresistance pxx increases linearly with B. 
From these experimental results, it is possible to plot the 
ratio LOT and consequently to extract the scattering rate 
1/r (figure H])). This clearly indicates that the scatter- 
ing rate linearly increases with the magnetic field. We 
note that there is also experimental evidence from opti- 
cal measurements that the scattering rate increases with 
increasin g rn agnetic field in these systems as it does in 
graphite [l3l . Since lo is given by huj — ^^^, and is so 
linear with B, we can deduce from experiments: 
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where 7 is a positive constant. Note that 1/t{B) > 
1/t{B = 0), and ljt < I/7 (with 7 equal to 0.4, 0.45 



and 1.25 for the present lfj,m x 5p,m, 0.27p,m x 6fj,m and 
100pm X 1000/iTO samples). Moreover, we clearly see (fig- 
ure HI)) that ujT is directly related to the system's ability 
to exhibit Shubnikov de Haas oscillations: At low ujt, the 
system does not exhibit such oscillations, while higher 
values are concomitant with the appearance of well sep- 
arated Landau-level in the electronic structure. A scat- 
tering rate l/r linear with B can thus explain both the 
unusual magnetoresistance behaviour and the quenching 
of Shubnikov de Haas oscillations. If one assumes that 
the coupling between the different planes is negligible it 
might seem quite difficult to explain a scattering rate 
that varies linearly with the magnetic field. Indeed, the 
conventional theory only predicts a quadratic correction 
to the magnetoresistance [18J. Yet, a scattering mech- 
anism due to coupling between the doped and undoped 
planes will depend on the magnetic field since the elec- 
tronic structure of the undoped planes depends on the 
magnetic field. Even in weak magnetic fields the Fermi 
level in the nearly neutral planes will be in the low in- 
dexes Landau levels and it will be in the zeroth Lan- 
dau level as soon as i? > Be = hnd/2e where Ud is the 
electron density. Infrared Landau spectroscopy for the 
nearly undoped planes indicate n^ is less than lO^^cm"^, 
i.e. Be « IT [11|. We will see that this scattering 
mechanism implies an increasing scattering rate with an 
order of magnitude of A{h/T{B)) ~ huj, with 7 close to 
1, compatible with the experiments. 

Model of electronic structure 

We consider that the perfect rotationally stacked 
planes are essentially decoupled since the coupling be- 
tween different states close to the Dirac point is of the 
order of ImeV (compared to 0.2 — O.SeV for A-B stack- 
ing) [12]. It is important to note that the hopping matrix 
elements between neighboring orbitals of the two rotated 
planes are of the order of 0.2 — 0.3eV, and that the effec- 
tive electronic decoupling is due to an averaging specific 
to the perfect structure. This averaging is destroyed by 
disorder, and consequently the defects in one plane will 
introduce a scattering in that plane and a coupling with 
electronic states in the other plane. Defects that spread 
on both planes also introduce in-plane scattering and in- 
terplane coupling. Here we consider a model with one 
doped plane and one plane which is essentially neutral, 
for which the zeroth Landau level is half-filled at all val- 
ues of the magnetic field. The two planes are decoupled 
in the absence of disorder but are coupled by disorder. 
In the following plane" 1" is the doped plane and plane 
"2" the undoped one. 

We consider also that the electronic states in plane " 1" 
are coupled essentially to the zeroth order Landau level of 
plane "2" only. Indeed we note that for B = IT the index 
of the Landau levels of the doped plane is n ~ 30 — 40 
and there are about 2Nl — 15 Landau levels of the doped 
plane which are closer to the zeroth Landau level than 
to the other Landau levels of plane " 2" . In the following 



we shall treat N^ as a large number. 

The model for the Green's functions of the uncoupled 
and perfect planes "1" and "2" are Gifi[z) and G2,o(-z) '■ 
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Fermi wavevector of the doped plane and uq its density of 
states at the Femi energy, without magnetic field. In this 
model we assume that the Landau levels of plane " 1" are 
equally spaced by Tilj where lo is the cyclotron frequency 
(which is valid as long as their index is high) . ELt) is the 
energy of the zeroth Landau level (i?Lo << NlTilu varies 
with B so that the zeroth Landau level stays half- filled) , 
and R is the ratio between the number of states in the 
zeroth Landau level of plane " 2" and in a Landau level of 
plane " 1" . R is equal to 1 if the two planes are equivalent 
(apart from the doping). We can also simulate the case 
where the doped plane is coupled to two undoped planes 
by taking R — 2. The degeneracy of a zeroth Landau 
level of a bilayer in AB stacking is also R = 2 [19| . Thus 
different configurations can lead to different values of R 
and we let this parameter free in the following. 

In order to treat the effect of in-plane scattering by 
disorder and intra-plane coupling by disorder we use a 
standard Self-Consistent Born Approximation (SCBA). 
We introduce the Green's function Gp{z) and the self- 
energies T,p{z) of plane "P" (P = 1 or P == 2). The 
density of states per unit surface in plane P is given by 
np{E) = —l/'Klm{Gp{z = E + ie)), where e is an in- 
finitely small positive real number. One get: 



Gp(z) = Gp,o(z-I]p(z)) 



Ep(z) = |ypGp.(z)-i 



2rp, 



(3) 
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P' is the plane coupled to plane P (P' = 1 if P = 2 
and vice- versa) . V'^ is an average value of the square cou- 
pling between states in plane "1" and in plane "2". The 
terms h/2Tpp represent the effect of in-plane scattering 
for each plane P (rp^p is the in-plane scattering time ). 

In dH El HI) the two important parameters for the effect 
of the interplane coupling are R and V. We analyse below 
the regimes of large and small R or V. We show that in 
the large V limit the experimental magnetic field depen- 
dence of the scattering rate (see figure [T]) are explained 
with a reasonnable value of the parameter R . We em- 
phasize that the physics described below presents some 
analogies with the coupling between localized d-orbitals 
with extended sp-states [20|. 

Large and small R regimes 



In this part, we focus on the case z = Elq + ie. Since 
the zeroth Landau level is half-filled, Elq is always close 
to the Fermi level. Let us assume that Gi(z — Ej^q + 
ie) ~ —innQ and thus ni(£'io) — '^o ("^6 show below 
that this corresponds to the large R limit). After ^ 
[21 15), this occurs for \ImY,i{ELo + ie)\/huj » 1. In 
this Hmit, the real parts ReG2{ELo + ie), RcGi^Elo + 
ie), ReT,2{ELQ -\- ie), Re'Ei{ELo + ie) are all neghgible. 
Using the correspondence h/r — 2ImY,, where r is an 
electron lifetime, one may write the SCBA equations in 
a form similar to the Fermi Golden Rule i.e. — = — — \- 

2TTV^np> with ni ~ no and n2 — RN{B)^ the densities 
of states at z = Elq + ie and then: 
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Here fi/T2.2 and 2TTV'^nQ are respec- 
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tively the width of the zeroth Landau level due to disor- 
der in plane 2 and to coupling with plane " 1" in the limit 
where its density is hq. One sees that if disorder in plane 
"2" (term h/T2^2) increases then the scattering rate Ti/ti 
decreases. Indeed the scattering by plane "2" is favored 
by a strong density of states in plane " 2" , whereas the 
term fi/T2.2 tends to decrease this density. If a >> 1, the 
coupling between planes " 1" and " 2" has essentially no 
effect (i.e. loti ~ ^^ti,!) but in the opposite limit a << 1 
the scattering rate for electrons in plane "1" increases 
linearly with the magnetic field and luti < Tr/2R. 

We show now that one get the regime Gi(z = Elo + 
ie) ~ —inno at sufficiently large values of R. For 
simplicity we consider the limit ?i/t2,2 — ^/''"i,i — 0; 
but let us just note that h/ri^i and Ti/t2^2 have op- 
posite effects, since Ti/ti,i favors the large R regime 
while ?)./t2,2 favors a small R regime. The SCBA 
equations can be written with dimensionless quantities 
FP = M±;,^t^,'-^^,Rl where Pp = ^^£1 or 
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This equation is 



independent of the coupling parameter V, and in that 
case ni/no and 2Im'E,i/hLu = \/liJTi are functions only 
of Elq/Tioj and R. 

We consider ELo/fuj — and Elq/Tiw = 1/2 for which 
by symmetry Elq = Ep and the real parts of the self- 
energies and Green's functions are zero. The result is 
shown in figure ^. For large R typically R > 1.5 — 2 
one has ni{Ep)/no ~ 1 which is the criterion for the 
strong scattering regime. In that case one recovers the 
strong scattering limit H/ti ~ 2R/'!Thuj that is loti — 
tt/2R. Note also that in the large R regime the results are 
essentially independent of Elo/Tiuj. At R — 1 {i.e. not 
in the large R regime) and for Epo/Uuj = the density 
of states ni{Ep) diverges. Indeed as soon as P < 1 there 
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FIG. 2: Left: value of ni{EF)/no as a function of R for 
ELo/fiuj — (full line) and for ELo/hu> — 1/2 (dashed line). 
Right: Same for the ratio cjt. 
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FIG. 3: Left: Dimensionless density of states in plane "1" 
{ni{E)/no) as a function of energy for different coupling 
V'^no/huj = 1 (dashed dotted line); 6.10"^ (fuU line); 1.10"^ 
(dotted line) and R = 2, Elo/Tilo = 1/2 . Right: same for the 
dimensionless density of states {n2{E)/no) in plane "2". 



are less states in the zeroth Landau level of plane "2" 
than in Landau levels of plane " 1" . This means that 
there exist uncoupled states of the Landau level of plane 
"1" and thus an infinite density ni(i?i?) as soon as i? < 1. 

Large and small V regimes 

We define the large and small V regimes respectively 
by 2'KV'^no/hu; » 1 and 2TrV'^no/hLj « 1. After the 
dimensional analysis the width W of the zeroth Landau 
level of plane " 2" and thus the energy range on which the 
electronic structure is modified by the coupling satisfies 
W/hto = W{ELo/huj, V^no/hiv, R). 

In the large V regime, as long as |z — ELo\/fiu! < 
2TTV^no/hiL; the term \z — Elo] « \V^Gi{z)\ and one 

The dimension- 



N{B) 



recovers Gi{z) — Gifl{z — R 

less Gi{z)/nQ and Yii{z)/fiu) depend on z/Ulo and R 
but not on Elq/Tlu) and 2ttV'^ nQ/hio. Note that the 
periodicity Tilo of Gi^q{z = E + ie) ([2]) implies the 
same periodicity for Gi{z — E + ie) and ni(E) in this 
limit. For the density of states in plane "2" one has af- 
ter ja El [l n2{E)/no = iR/TT^)ihuj/V^na)ino/ni{E)). 
Thus n2{E)/no presents the same periodicity (within 
the range W) and is very small. The spectral weight 
RN{B) of the zeroth Landau level spreads on a width 
W and thus since n2{E)/no ~ (i?/7r^)(?iw/F^no) one 
has W/hui ~ Tv'^V'^no/hLU. Finally in the smallF regime 
2TTV'^no/huj << 1 the width W/hco << 1 and the Landau 
levels of plane "1" and "2" are only slightly hybridized 



(see Figure ^). 

Magnetotransport 

When the density of states is uniform on an energy 
scale W >> hui, which is the case in the large R and 
V regime (see figure ([3])), we can apply the semiclassical 
theory of transport. The scattering time ti is given by 
H/ti — Ti/ti^i + 2hujR/{n{l + a)) where Tj/n^i is the 
in-plane scattering rate. As long as a is not too large 
we expect the term 2hu)R/{'K{l + a)) to be of order Tiuj 
and the model is consistent with the experimental results 
presented above (see figure[l]and equation[T]). Shubnikov 
de Haas oscillations can occur when the field dependent 
scattering studied here is destroyed, that is for a >> 1. 
This can be due for example to disorder in plane "2" or 
even to a confinement effect as in a ribbon of finite width 
|2lj . as in figure ([1]). Shubnikov de Haas oscillations can 
thus be enhanced by disorder or by confinement effects. 
This spectacular effect is clearly observed in experiments 
4]. 

To conclude we have shown that magnetotransport ex- 
periments on epitaxial graphene are consistent with a 
scattering time that is magnetic field dependent and is 
reduced to the order of the cyclotron period. This ex- 
plains the unusual variation of the magnetoresistance, the 
quenching of the Shubnikov de Haas oscillations and of 
the quantum Hall effect . The magnitude of the field de- 
pendent scattering time is consistent with a mechanism 
where the conducting electrons of the doped plane are 
scattered due to their coupling with the zeroth Landau 
level of the undoped planes . 
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